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Discrete analogues in harmonic analysis: 
Spherical averages 

By A. Magyar, E. M. Stein, and S. Wainger* 
Abstract 

In this paper we prove an analogue in the discrete setting of of the 
spherical maximal theorem for M. d . The methods used are two-fold: the appli- 
cation of certain "sampling" techniques, and ideas arising in the study of the 
number of representations of an integer as a sum of d squares, in particular, the 
"circle method" . The results we obtained are by necessity limited to d > 5, and 
moreover the range of p for the IP estimates differs from its analogue in M. d . 



1. Introduction 

Geometric considerations, in particular curvature, play an important role 
in harmonic analysis in M. d . Emblematic of this are the properties of the 
spherical maximal function. Given the significance of this operator, it is an 
interesting and natural question to ask what happens when we consider its 
discrete analogue; that is, what can be said of the corresponding version of the 
spherical maximal theorem taken over It is the purpose of this paper to 
answer this question by proving optimal £ p estimates in this setting. 

We shall now describe these results, turning first to M d . The spherical 
averages are defined by the operators A\,wheve 

(/) = /* da x 

with da\ the normalized invariant measure on the sphere |x| = A. With the 
definition of the maximal function, A*(f)(x) = sup \A\(f)(x)\, we recall the 

main estimate for it, 



0<A<oo 



1.1) || A* (/) < A ||/|| LP(Rd) , if P>-^zj and d - 2 ' 



(See [S], [SW1], [Bl].) 



*This work was supported by NSF Grants: DMS-9970899, DMS-9706889, and DMS-9731647. 
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The discrete analogue of A\ is the operator 
(1-2) A x (f)(n) = j±- £ f(n-m). 

|m|=A 

Here n and m are restricted to range over Z d ; also N(X) = the number of 
m € 1j d , so that \m\ = A. Notice that only those A for which A 2 is an integer 
are relevant; also observe that N(X) = r^(A 2 ), where r^(/c) is the standard 
counting function giving the number of ways of representing A; as a sum of d 
squares. 

Now, up to this point, formulating a discrete analogue of the spherical 
maximal function, i.e. A*(/)(n) = sup \A\(f)(n)\, and asking the question 

0<A<oo 

of its £ p boundedness, seem quite straightforward.* However, this is misleading 
since quite different ideas must come into play in the discrete analogue, and 
anyway, the range of exponents is not the same as the version in R d . The 
theorem we prove is the following optimal result. 

Theorem. The maximal operator A* is bounded in t v (JL d ) to itself for 
p > > when d > 5. 

Alex Ionescu has pointed out to us that simple examples show that this 
result cannot be improved: in fact, when d > 5, A* is not bounded on t v for 
p < moreover when d < 5, the t v boundedness fails for every p < oo 

(the case p = oo is of course trivial). The relevant examples can be found 
in Section 8. Here the facts that the number of representations ra(k) is an 
irregular function of k when d < 4, while r^(/c) ps k~^~ when d > 5, play a 
role. (For these assertions about rd(k), consult [W].) 

Our attack on the discrete spherical maximal function proceeds in three 
stages. To begin with (motivated by the ideas of the circle method) we ap- 
proximate A\ by an infinite sum of simpler operators 

(1.3) M x = c d Y< e- 27Tix2a / q M a x /q , 

with each M^ q associated to a reduced fraction a/q, with < a/q < 1. Now 
since each M^ q is a convolution operator on Z d , it corresponds to a Fourier 
multiplier m^ / '' ? (£), which is given by 

< /q (0 = E G(a/q,l)<t> q (£ -i/q) da x (Z -i/q). 



* Here we use the notation that for / defined on Z d , it belongs to £P(Z d ) if \J |/(n)[ p is finite. 
The i v norm is of course the p th root of the sum. 
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Here G is a normalized Gauss sum, § q is a suitable cut-off function, and da\ 
is the Fourier transform of the unit measure da\ on the sphere \x\ = A. 

Notice that the first term of M\, corresponding to a/q = 1 = Omod 1, can 
be viewed as the vestige of the continuous analogue on M. d . All the other terms 
are approximations corresponding to the other rationals. 

The second stage is to study each M^ q as a sort of discrete analogue of 
an operator on R d . The main tool is a general abstract theorem which allows 
one to pass from certain convolution operators on M. d to analogous operators 
on 7h d . While ideas about special cases of this principle have been implicit in 
the past, our general approach seems both new interesting in its own right. 
It is presented in Section 2. It is based in part on variants of "sampling" 
ideas which go back to Plancherel and Polya [PP] and which were taken up 
again later by Shannon [ShW]. Using arguments of a different kind, Bourgain 
obtained certain results of this form; see [B2, (3.5)]. 

The final stage of the argument is to show that M\ is an adequate approx- 
imation of A\. This is begun in Sections 4 and 5, and is concluded in Sections 
6 and 7. 

The analysis of our theorem has as its starting point a partial result ob- 
tained previously by one of us [M] (see Proposition (4.2) below). The interested 
reader may also want to compare the related ways the sums V] g 2™-£ are 

\n\=\ 

treated in our paper (see Section 5), and in a previous work of Bleher and 
Bourgain [BB, §6]. The context of that paper is however quite different from 
ours. 



2. Discrete analogues of convolution operators 

Suppose T(f) = f * K is a convolution operator in L p (M. d ) to itself with 
a suitable distribution kernel K. Then, as is known, its Fourier transform 

K = m(0 = f K( X )e-^d X is a bounded function, and we can think of T 

jR d 

as a Fourier multiplier operator given by (Tf) A (^) = m(^)f(^). 

To be precise, in what follows we shall assume in this section that in 
addition to m(£) being bounded, it is supported in the fundamental cube 
Q = {£_ = : -1/2 < £j < 1/2, j = !,...£}. In this case K(x) = 

/ e 27m ^ m(ti)dti is an L 2 function on M. d , which is continuous (in fact, C°°). 

Thus Kdis = K is well-defined, as is the convolution operator acting on 

z 

functions on Z d given by 

2dis(/) = Z*ifdis, T dis (f)(n) = Yl K(m)f(n-m). 
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Note that the condition that the multiplier be supported in Q is natural. 
Because then not only does T determine T^ S1 but conversely T^ s determines 



T, i.e. K 



determines K. This follows since K 



determines the Fourier 



coefficients of the function when expanded as a function on Q. 

Let m per be the periodic extension of m, i.e. m per (^) = E m(£ — £). 



Then m per (£) is the Fourier multiplier corresponding to in the sense 
E T dis (f)(n)e- 2m < = m pcr (£) ]T /(n)e" 2 ™«, 



for suitable functions / on Z d . 
Let us note that 



(2.0) m per (£) = £ tf( 



in the sense of L 2 convergence of the series on any compact subset of R d . 

In fact, m(£) = V] K(n) e ~ 27rm £ on q represents the Fourier inversion 

of the identity K(n) = J m(^)e 2nm ^ d£, (m(£) is supported in Q); and, more- 
over, m per (£) is the periodic function which agrees with m(£) on Q. This 
establishes (2.0). 

The question we will be concerned with is how the norm of T dis as an 
operator on £ p (Z d ) is controlled by the norm of the operator T acting on 
L p (R d ). For our applications it will be important to be able to deal with the 
more general case where the L p and t v spaces of complex- valued functions are 
replaced by the spaces L^(M. d ) and £ p B CZ d ) of functions taking their values in 
the Banach space B. In order to avoid technical problems involving definability, 
measurability, etc., we shall restrict our attention to the case when the Banach 
spaces in question are finite-dimensional. However, all our estimates will be 
independent of the Banach spaces in question, so that a limiting argument will 
encompass the results in the generality needed. In particular, this argument 
will apply to the case when B is an L°° space, which is what is needed for the 
maximal theorems below. 

We shall suppose that B\ and B2 are a pair of finite-dimensional Banach 
spaces, and assume that m(£) is a bounded measurable function, taking its 
values in C(B\, B2); and as we have said we suppose m is supported in Q. 
Then T, described above, is a bounded mapping from L 2 Bl (R d ) to L 2 B2 (R d ), 
and similarly T dis is bounded from i\ (Z d ) to £ 2 B2 {% d ). 

PROPOSITION 2.1. Fix p, 1 < p < 00. IfTis bounded from L p Bi (R d ) to 
L p B2 (R d ), then T dis is bounded from e p Bl (Z d ) to £ p B2 (Z d ). For these operators 
we have the norm inequality 



DISCRETE ANALOGUES IN HARMONIC ANALYSIS 



193 



(2.1) T dis < C T 

v ' uls pP .pP — TP . rP ' 

l B 1 ~* l B 2 lj B l ^B 2 

with a bound C that depends only on the dimension d, but not on p or the 
Banach spaces B\ and B2. 

Remarks. (1) It would be interesting to know if C can be taken to be 
independent of the dimension d, or for that matter if C = 1. 

(2) There is a converse to (2.1); i.e., a reverse inequality also holds. Since 
that fact will not be used below we will omit its proof. 

The proof of the proposition requires the following "sampling" and exten- 
sion lemma. We fix the function \& on M. d by 

/sin7rxi\ 2 /sin7TX2\ 2 f smTTX d \ 2 
tf(ar) = ... , x = (xi,x 2 , • • .x d ) . 

\ TTXi J \ 1TX2 J V TTX d J 

For any suitable function / on Z d we consider its extension / ext = F on R d 
given by 

(2.2) F(x) = / cx t(x) = Yl f(n)V(x-n). 

(Note that if / € £ p for some p, the series above converges for every x G M. d .) 
We observe that in fact F = f, since ^(0) = 1, and ^(n) = 0, if n G Z d , 
n ^ 0; thus f cxt is a genuine extension of /. The following estimate holds for 
any (finite-dimensional) Banach space B. 

Lemma 2.1. // / € F(Z d , B), then F e 17 (R d , B), and 

(2-3) (1M) 11/11,, < ||F|| l p < A 11/11,, . 

Here A is a constant that depends only on d, but not p or the space B. 

Ideas of this kind go back to Plancherel and Polya [PP]. In that work 
(when e.g. d = 1), the function S1 ^ x was used in effect in place of ( S1 "^ a ) • 
The resulting version of (2.3) is then more delicate and holds only in the range 
1 < p < 00, since it involves the Hilbert transform; it also does not cover the 
case of Banach space- valued functions. 

To prove the lemma we observe two easily established estimates, 

/ ^(x)dx < Ax and sup V ty(x - n) < A 1 . 

jRd X 

Then for any p < 00, by Holder' inequality 

|/ext(x)| p < \f(n)\*9(x-n)) fe*(x-n)V 
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and integration in x then gives 

II f + Il p < II f ll p A p 
|| Jext \\L V B — \\ J \\£ B 1 ' 

The proof of the corresponding result for p = oo is similar but simpler. 

To prove the converse inequality, choose $ to be a C°° function with 

compact support so that $(£) = 1, when £ € 2Q. Since 

(1 — |£l|)+, it follows that 

= and hence VP * $ = * . 
As a result / cxt * $ = / cxt and since / ex t(n) =/("), we have 



/ sin 7rxi \ ' 



f(n) 



/ext(y) < &(n-y)dy. 



Thus as before, 



< 



$(n-y) 



and 



if 



/ext(y) $(n - y) 

ll/II^P < ll/cxt|| L P ^2' 

/ \$(y)\dy < A 2 , and sup V |$(n - < A 2 . 



p-i 



The argument also gives the case p = oo. The lemma, inequality (2.3), is thus 
established with A = max(Ai, A2). 

To prove the proposition, we consider the cube 3Q which can be cov- 
ered by 3 d disjoint translates of Q. In fact, it is easily verified that 3Q = 
U (Q +£). Now let m(£) be continued periodically to 3Q, i.e. de- 

BUp|^|<l 

j 

fine m(£) by m(£) = V] m(£ + £). We let T denote Fourier multiplier 

SUPj|^j|<l 

operator, whose multiplier is m(£). Then clearly 

3 d 



(2.4) 



< 3° 



B 2 



On the other hand, we claim that 
(2-5) f (/ cxt ) = (T dis (/)) ext . 

To verify (2.5) it suffices to do it for / = 5 m , for every fixed m, where 

1 if n = m 



$ m (n) 



if n 7^ m . 
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We will check this by taking the Fourier transform of both sides of (2.5). 
Indeed 

T dis (f)(n) = K(n-m) 



and 



Hence, 



(?dis(/)ext) / 



(T dis (/))ext = E K(n-m)*(x-n) 



E K ( n - m) § (0 e~ 2 ™ ? 

n 

(E K(n - m) e ~ 2 ™«) 

f E * (») e 



= m pcr (0^(0e- 2mm « (by (2.0)). 

On the other hand 

T(/ext) A = m(0(/cxt) A = m(0(^(x-m)) A = S(#(0e~ 2 ™«. 

Now we have the desired identity, since tti(^) — Tfiper (£) on the support of ^ 
(note that 2Q C 3Q). 

Once (2.5) is established we have 



?di s (/) 



< A 



B 2 



(7dis(/)) 



ext 





n/ext) 


< 


r(/ext) 











< 3 d ,4 

< 3 d ,4 2 



T 
T 



L Sl^ L S 2 



/ext 
/ 



(by the lemma) 

(by (2.4))) 

(by the lemma). 



Thus 



C = 3 d A 



Tdv. 



d A 2 



< 3 d A 2 



pP .pP 



, and the proposition is proved with 



We now fix an integer q > 1. We shall also make the stronger assumption 
that m(£) is supported in Q/q, and consider m^ cr defined by 



(2.6) 



mlM = E m^-l/q). 
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Notice that m q ev is periodic with respect to elements in (l/g)Z d , and 
hence, in particular, periodic with respect to Z d . 

We consider the operator Tj is , a convolution operator on Z d , having m q er 
as its Fourier multiplier; i.e., 

E ^ is (/)(n)e- 2 ^ = ml er (0 £ /(n) e - 2 ™-5, 



for suitable /. 

Corollary 2.1. 

(2.7) 



J dis 


< C 


T 









Again the bound C does not depend on p, B\ and B2; it is also independent 
ofq- 

Proof of the corollary. Let T q be the operator on L P B (W 1 ) to L B 
whose multiplier is m(£/q). Notice that m(£/q) is supported for £ G Q. 
Now a simple scaling argument shows 



rpq 



L Si^ L s 2 



L Si^ L s 2 



and so if (T q )^ is is the discrete analogue in the sense of Proposition 1, 



(2.8) 



(T q ) 



(lis 



< C 



B 1 l B 2 



However, we must emphasize that (T q )di s 7^ T3- . In fact, the convolu- 
tion kernel of (T 9 )^, which comes from the multiplier m(£/q), is K q (n) = 
q d K(qn), n G Z d . 

Next we observe the convolution kernel, K#(n), of T^ is is given by 



K#(m) = < 



q d K{m), ifmeqZ d 

= if m G Z d , but m <£ qZ d , 



because 



= J ( E e 2 ^ m de 



7 d if m G (/Z d 



if m $ qL d . 
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Now finally let T# denote the operator mapping functions of qL d to itself, given 
by the kernel K#, i.e. 



T#(f')(nq) = E f'((n-m)q)K*(mq) 



Then clearly 



'-# 



>4 2 («z d ) 



dis 



-B 2 



gZ , given by 



which is an immediate consequence of the isomorphism 7L a 

n «-> gn, (n G 

Finally note that T% is can be written as T # ® J, if we write ^(Z d ) 
as l p B (qL d ) ® l v {jL d /q7L d \ with T# acting on the first factor, and the iden- 
tity acting on the second factor. 

As a result 



J dis 


< 






l B 1 ^ l B 2 





B 2 

Combining this with (2.8) proves Corollary 2.1. 

We next consider a version of a convolution operator, whose multiplier is 
somewhat akin to (2.4). Here we shall consider 

(2.9) m(0 = 

£ez d 

under the following assumptions: 

(a) is a C°° function supported on Q/q. As a function on Q it has the 

Fourier expansion 

2nimt; 



with 



m = E <p« 

E \frn\ < A. 



(b) {7^} is a <7^ d periodic sequence; i.e., 7£ = 7^ if ^ — £' G (?Z d . 
Now let {-%} be the Fourier transform of {7^}; i.e., 7 S = 



e 2jnst/q Je 



E 

We shall also restrict our attention to scalar-related functions on £ p (7j d ), as 
opposed to the Banach-space case treated in the previous proposition, because 
of the specific use of PlancherePs identity. Our result is as follows. 
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Proposition 2.2. Let T be the operator on functions on Z d whose 
Fourier multiplier is given by (2.8), satisfying the conditions above. Then with 
1<P<2, 

/ \2-2/p / \2/p-l 

(2.10) T < A sup fa | sup \% | 

£p(Z d )^£P(Z d ) W / V * / 

For the case p = 2 we use Parseval-Plancherel's theorem for together 

with the disjointedness of the supports of the $(£ — i/q) and the fact that 

sup|<I>(^)| < A. This implies that |m(£)| < Asup|7^|, yielding the case p = 2. 
C ' < 

For the case p = 1, we calculate the ^(Z d ) norm of the kernel K(n), 

corresponding to the multiplier m(£). It is given by 



J Q 



Hence by property (a), 1-^(^)1 < A sup \j n |, and as a result the case 

p = 1 of (2.10) is proved. The general result for 1 < p < 2 then follows by 
Riesz' convexity theorem. 

3. The main term 

The averages we are interested in, 

Mf)(n) = j^r £ /(«-"»), 

|m|=A 



will be replaced by the equivalent averages 

1 



J2 f{n-m) 



\d-2 

H=A 

when d > 5. This equivalence comes about because, as we have pointed out 
with N(X) = number of n € so that |n| = A, we have N(X) ~ \ d ~ 2 , 
whenever A 2 is an integer, and d > 5. In order not to introduce new notation, 
we shall designate these averages also by A\ and now write 

(3-1) A x (f)(n) = -1_ ]T f(n-m), 

\rn\=X 

and in what follows we shall always assume that A is restricted so that A 2 is 
an integer. 
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Here we shall deal with the main term in the approximation of A\. It is 
a convolution operator M\ acting on functions on Z d , which can be written as 

oo 

(3.2) M x = c d Y: E e~^ x2a ^M a x /q , 

q=l \<a<q 
(o, 9 ) = l 

where the sum is taken over all reduced fractions a/q, with < a/q < 1. Hence 
Cd is the constant = y(1/2) ■ Also each M^ q is the convolution operator whose 
multiplier is 

(3.3) E G (%> * 9 (£ - V?) ddx (£ - */<?) • 



In the above, ^g(C) = ^(<z£)> w here ^ is a (7°° cut-off function supported 
in the cube Q/2, with *(£) = 1, for £ € Q/4. Also G(a/q,£) is the normalized 
Gauss sum 

d e 27ri(|n| 2 a/q+n^/q) 



(3.4) G(a/ g ,£) = g- d E 



and d&\(£) is the Fourier transform of the normalized invariant measure da\ 
supported on the sphere \x\ = A. Note that (3.3) is periodic on £ with periods 
in Z d , since G(a/q,£) = G(a/q,£'), if £ = £ f modqZ d ; also, for each £ only one 
term in (3.3) is nonzero. 

We define the corresponding maximal operators, 



M*(/)(n) = sup |M A (/)(n)|, 

0<A<oo 



and 



M a J q (f){n) = sup M^(f){n) 
0<A<oo 

The basic estimates for these are as follows: 

Proposition 3.1. (a) || Ml lq ||^_^>= O (g-^ 1 - 1 /?)) if d > 3, and 
(b) || M* ||#>-#> < A if d > 5, and ^ < P < 2 • 

To prove part (a) we write ^ = ^ • where \P' is another C°° function, 
supported in Q, with = 1 for £ G Q/2. Then the operator corresponding 

to the multiplier (3.3) can be written as a product of two operators, with 
multipliers respectively: 

E G(a/q,e)%(£-e/q) 

and 

E * q (S-t/q)d*x(S-e/q), 



T a/q, 
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where = *&'(qO> ^ we reca U that for each £ only one term in each of the 

above sums is nonvanishing. 

To the first multiplier we apply Proposition 2.2 (in §2) with 7^ = G(a/q, £), 
and $>(£) = Note that $(£) is supported in Q/q, and moreover 

¥>(n) = / <&(0e- 2 ™^ = g- 1 *^), 

where \& is the Fourier transform of Now \*ff(x) \ < An (I + |#|) _iV for all 
N > 0, so that ^ |p(n)| < A. 
nei d 

Next, there is the estimate \G(a/q,l)\ = 0(q~ d / 2 ); this is well-known, 
but in any case it follows from the standard one-dimensional case merely by 
observation that G(a/q,£) is a d-fold product of these one-dimensional sums.* 
Moreover, if % = £ ^is-e/q G ( a / q ^ ^ then> 



Af g = J_ V e 2^^/<i, e 27rj(a/<i,)|n| 2 ^-Kin-ijq = ^{a/ q)\s\ 2 _ 

9 nez d / qZd £ez d / q z d 

Hence by Proposition 2.2, the norm of the corresponding operator (acting on 
F to itself, 1 < p < 2, with scalar- valued functions) is 0(g-( d / 2 )( 2 - 2 /p)) = 
( 9 -d(i-i/ P ))_ 

Next, the multiplier 
(3.5) £<Mf-*/9)d<7A(£-V?) 



corresponds to a convolution operator from £ p (Z d ) (scalar- valued) , to £ p B (Ij d ), 
where B is the space of functions of A > 0, for which A 2 is an integer, and 
< A 2 < N. Notice that = is a bounded multiplier of L p {R d ) 

to itself (with norm independent of q). Observe also that dd\{£) is a bounded 
multiplier from L p (R d ) to L B (R d ), for p > ^4r> w hich is a consequence of the 
spherical maximal theorem in M, d . Finally, note that m(£) = & q (£,) d&\(£) is 
supported in Q/q. Thus, applying the corollary to Proposition 2.1, we see that 
(3.5) is a bounded multiplier from £ P (/L d ) to £ B (l* d ), with norm independent of 
-/V (and q). Letting N — ► 00, and combining this with the estimate for the first 
multiplier, we have established conclusion (a) of Proposition 3.1. The second 
conclusion follows from this because 



m* <c d J2 J2 M - 



ajq 



\<q<CO (a,g)=l 
l<a<q 



tFor the one-dimensional estimates, see [W]. 
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so that 

< oo , 

^ p 1< 9 <oo 



if 1 - d(l -1/p) < -1, i.e. when p > 



4. Approximations 

We now state the assertions which guarantee that M x provides an ade- 
quate approximation to our operator A x . There are two facts; the first is a 
purely £ 2 statement. 

PROPOSITION 4.1. There is a bound A, so that for any A > 0, 

(4.1) || sup \A x (f) -M x (f)\ y< AA 2 ~ d / 2 \\ f\\p, ifd>5. 

A<A<2A 

The second is a partial result for A\ which was known previously (see [M] ) . 
Proposition 4.2. There is a bound A, so that for any A > 

(4.2) || sup \A x (f) We, < A || / \\ £P , if d > 5 , p > -J— . 

A<A<2A a - I 

Recall that the A which appear in (4.1) and (4.2) are always restricted to 
the fact that A 2 is an integer. 

We shall momentarily take these two propositions for granted and see how 
they, together with Proposition (3.1), prove our main theorem. 

Now (4.2) together with Proposition (3.1) yield 

|| sup (A x - M x )f Hp < A || / || £P for 2 > p > . 

A<A<2A d — I 

Interpolating this with (4.1) gives 

|| sup \(A X -M x )(f)\ \\ lP < AA-z(p) || / ||^ 

A<A<2A 

for some s(p) > 0, if < p < 2. 
Next, 

oo 

sup \(A X -M x )(f)\ < sup \(A X -M x )f\. 

l<A<oo k=Q 2 k <X<2 k + 1 

Taking the t v norm we get that 

|| sup \(A X - M x ) f\ Wp <A'\\f , for -J— < p < 2, 
KA<oo a — 2 
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since E^ ^ p ^ k < oo. Thus, invoking Proposition (3.1) again yields 

k 

II sup \A x (f)\ \Up<A || / \\ip when — ^— < p < 2. 

i<a<oo a — I 

Since the corresponding estimate for p = oo is trivial, the full range ^2 < 
p < 00 then follows by interpolation, proving the main theorem. 



5. The decomposition of A\ 

To prove the crucial approximation property (4.1) we shall decompose 
the operator A\ into a sum, each of whose terms corresponds to a fraction a/q, 
with 1 < q, 1 < a < q, and (a, q) = 1. It is here we use the ideas of the "circle 
method" of Hardy, Littlewood, and Ramanujan. 

Let us fix A > 0, and consider any A for which A < A < 2A. We shall 
write oa(£) for the multiplier corresponding to the operator A\ given by (3.1). 
We claim that 



2tts\ 



2 



A . Jo 



Here e is positive, but otherwise arbitrary; we will fix it later by setting e = 

A 22 

1/A 2 . This identity is obvious because / e 2?r ^' n ' ~ A )* dt = 1 or according 

Jo 

to whether |n| = A or not. 

Now we introduce the function 



(5.1) F(z,0 = £ 



g— 2n\n\ 2 z g27rin£ 



for 3?(z) > 0, and we make a Farey direction of level = A of the interval 
[0,1] of the t integration. That is, for each a/q, (a,q) = 1 with 1 < a < q, 
and q < A, we associate the interval I(a/q) = ji : — < t — a/q < 
where a = a(a/q,A) sa 1, and (3 = (3(a/q,A) « 1, with a and /? chosen 
appropriately. We denote by 1 (a/q) the corresponding intervals translated to 
the origin, I(a/q) = |r: — ^ < r < p-j. Inserting this in the above formula 
for a\(£) and using identity (5.1) we get 



a x 

where 

a 2ire\ 2 



(0 = E E 



1<Q<A l<a<? 

(a,«)=l 



(5.2) a7 9 (0 = e - 2 ™ A V9 / 6 -2™a 2 t _ - T _ • , j ^ dr _ 

A" - J/(o/g) 
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Next we use the fundamental identity for the O function (5.1). It states 
that for 9ft(z) > 0, 

Here G(a/q, £) is the normalized Gauss sum (3.4). The above is the cf-dimensional 
version of a familiar identity. (For d = 1 see, e.g., [SWi, (3.4)]; also [W].) The 
general case d > 1 can be proved the same way invoking the Poisson summation 
formula; alternatively one can observe that (5.3) is merely the ci-fold product 
of the corresponding 1-dimensional identities for each variable £i, £2, • • • £<2j sep- 
arately. 

From (5.3) and (5.2) it follows that 

(5.4) a a x /q (0 = e -2«AV, £ G (a/g, £) J A (a/ 9) $ - i/q) , 

where 

(5.5) J A (a/ g ,C) = -T-fY- / e" 2 ™ A T (2(e - *T))" d / 2 eW^J dr . 

6. Approximations, continued 

We shall approximate the multipliers a x q {£ 1 ) above by multipliers b a x q {£) 
where the cut-off factors & q (£ — l/q) have been inserted in (5.4). That is, we 
define 

(6.1) b a x /q (0 = e" 2 ^^ G ( a /V> ') ~ Z/V) J A« " i/q) ■ 

Here = <%£)■ 

Next we approximate b a x q {^) by replacing the integral (5.5) that appears 
in (6.1) by the corresponding integration when taken over the whole real line. 
So we set 

(6.2) c a x /q (0 = £ G(o/(Zj £) ^ _ e/q) j a( ^ _ t/q) ) 
with 

271-eA 2 ^00 „ — 7t 1 § 1 2 

(6-3) h (0 = ^32- j /_ oo e- 2mX T (2(e - ir))^ 2 e 2 ^) dr . 

We define the operators ^f 9 , -B"^ 9 , C^ 9 , as the convolution operators 
(acting on functions of Z d ), whose Fourier multipliers are respectively, a°/ 9 (£), 
&f (£), and cf (£)• 
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Proposition 6.1. 

( 6 - 4 ) E E II su p 

Kg<A i<a.< q A<A<2A 

(a,q) = l 

( 6 - 5 ) E E II su p 

Kg<A i<«<9 A<A<2A 

(o,g)=l 



||^ <^4A 2 - d/2 ||/|| 

- c A a/9 ) /| y < Ah 2 ~ d / 2 || / 



It is understood that in the above assertions our e is fixed to be = 1/A 2 . 
To prove (6.4), let F T be the function on Z d which is given in terms of its 
Fourier expansion by F T (£) = where 

MO = E (! -^))e-^-^l 2 /2( £ —). 

Note that since each term in the sum is supported where |£ — £/q\ > c/q, 

—ce 



sup |/i(£)| < -4 exp 



for some c > . Thus 



g^ + r 2 ), 
F T \\ P < A exp ( g2(g2 ^ r2) ) II / \\p ■ 



Now observe that 



sup 

A<A<2A 



Jl(a/q) 



because G{a/q,£) = 0{q- d l 2 ). 
As a result, 



sup 

A<A<2A 



< AhT d+2 q- d ' 2 ( {e 2 + r 2 )- rf / 4 exp 

Jl(a/q) 



-ce 



q 2 {e 2 + T 2 ) 



dr- || / 



Now, because e u < cu rf / 4 , we get a contribution of A d+2 e d ^\I{a/q)\ 
| / Taking into account that e = 1/A 2 , and ^ ^ \I{a/q)\ = 1, 

1<<?<A a,q 

we obtain (6.4). 

The proof of (6.5) is similar. Notice that we are now integrating over r in 
the complement of I(a/q), and thus |r| > c/qA. We are led in the same way 
to see that 



sup 

A<A<2A 



(B*f* - C a x ,q ) f 

< AA~ d+2 q- d / 2 



[ r- d ' 2 dr- || / 

J\r\>c/qA 



< AK- d+2 q- d l 2 (qKf 2 - 1 || f\\p= AK- d ' 2+l q- 1 \\ f 
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Now sum in a, then over q, q < A. This gives a contribution of 

o(A- d/2+1 ) ( £ Ei) T 1 \\f\\e=o (A- d/2+2 ) II / y, 

\1<9<A l<a<g / 

which proves (6.5). 

To complete the approximation process (the proof of (4.1)) we now identify 
h(0 given by (6.3). 

Lemma 6.1. 

h(0 = c d da x (0. 

Taking this temporarily for granted we observe that as a result, c d M^ q = 
C a x ,q (see (3.3) and (6.2), (6.3)). Hence for A < A < 2A, 

\(A X -M x )f\ < E E IK' " C x q ) f 

1<9<A (a,q)=l 

+ E E *K 9 / 

<?>A (a,q)=l 



However, 



A a/q _ c a/q 



< 



At~B x 



a/q 



+ 



£>a/q _ (ja/q 
A A 



Thus, invoking Proposition (6.1), and Proposition (3.1) for p = 2, we see that 



sup | (A x - M\) (/) | \\p 

A<A<2A 



= o(a^) \\f\\ e2+ h £ 

\?>A l<a<q 

= o(A 2 - d / 2 ) ||/|U , if d > 5. 



-d/2 



fWt 



Therefore, Proposition (4.1) is now proved, and with it the proof of our 
main theorem is complete, save for verification of the lemma above. 



7. Proof of Lemma 6.1 

The identity 

p 2-rre\ 2 poo -t|«I 2 

(7.1) y_ Qo e- 2mA T (2(e - ir))^ 2 e^> dr = c d da x (0 

is probably known, but we have not found it in the literature, and so we will 
give a proof. 
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First we observe that the left-side of (7.1) is in fact independent of e, and so 
we may take e = 1/A 2 . We see that this follows by changing the contour while 
integrating the function F{z) = (2z)~ d / 2 e 2lTZ e ^ 1^1 l 2z along lines parallel to 
the x axis in the upper half-plane. Next, with e = 1/A 2 , and with the change 
of variables A 2 r = t, 

/OO 1 -irA 2 |£| 2 
e ~ 2nit 7 — 7 -—777 dt . 

(2(1 - it)) d / 2 

We now insert an extra convergence factor e _7r5 * 2 in the integral defining 
I\ above. Denoting the resulting integral by ijj we have ij[ — > Jj[; moreover if 
if is any test function in the Schwartz space, then 

fmhiOdt = iim f miKoda. 

Also, 

(7.2) [m4(0<% = [Mx)i s x (x)dx. 

JR d JR d 

— vr|A| 2 g| 2 

Calculating the Fourier transform of the Gaussian e 2 (»-«) we see that 
J — 00 

which in turn is e - 27r M 2 /A 2 ft- 1 / 2 e -n0--\ x \ 2 . Inserting this in (7.2), and 
letting 5 — > 0, we obtain 

/ h(0 d t = Cd <p(x)da x (x) , 

and thus = Cdda\(l;), as was to be proved. 

Note that 

d X{) loo (2(l-^)) rf / 2 T(d/2)- 



8. Counter-examples 

Since we shall be dealing with all d > 2, we return to the original definition 
of the averages A\, 

JV V A J | m | =A 

Let us take / to be the unit mass at the origin; i.e. /(0) = 1, and /(n) = 0, if 
n € Z d , n / 0. Then clearly / € £ p (Z d ), for every p. Next, we observe that 

= Vjv(a) = Vjv(H). if H = A. 
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Hence, 

(8.1) A(/)(n) = supA A (/)(n) > l/ N{H) . 

A 

(Recall that N{\) = number of m G so that \m\ = A ; i.e. 7V(A) = ni(A 2 ).) 

Consider now the situation when d > 5. Then as we have pointed out, 
N(X) w A rf ~ 2 , and so A*(/)(n) > c|n| _d+2 . But the latter function belongs to 
£ p (Z d ) only when p > and so the necessity of that condition is proved. 

Next assume d < 4. We shall use the fact that r4(2 2fe ) = 24, for every 
k > 1. This follows from the Jacobi formula which states that r 4(171) = 8 • 
cr*(m), where cr*(m) is the sum of the divisors of m which are not divisible by 
4. (See [HW, Chap. 20].) Thus 

r 2 (2 2fe ) < r 3 (2 2fc ) < r 4 (2 2fc ) = 24. 

Now for each d, d < 4, we then have N(X) < 24, if A = 2 k . And so for n € 7L d 
with n = (2 fc ,0, ...), we see that A*(/)(n) > 1/24, by (8.1). Because this 
happens for infinitely many n, we have A*(f) £ £ p , for any p < oo, and so the 
necessity of the condition d > 5 is established. 
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